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ABSTRACT 

We apply the ZTRACE algorithm to the optical NOG and infra-red PSCz galaxy 
catalogues to reconstruct the pattern of primordial fluctuations that have generated 
our local Universe. We check that the density fields traced by the two catalogues are 
well correlated, and consistent with a linear relation (either in 5 or in log(l + 5)) with 
relative bias (of NOG with respect to PSCz) 6rei = 1.1 ± 0.1. The relative bias relation 
is used to fill the optical zone of avoidance at \b\ < 20° using the PSCz galaxy density 
field. 

We perform extensive testing on simulated galaxy catalogues to optimize the re- 
construction. The quality of the reconstruction is predicted to be good at large scales, 
up to a limiting wavenumber fcii,„ ~ 0.4 h/Mpc beyond which all information is lost. 
We find that the improvement due to the denser sampling of the optical catalogue is 
compensated by the uncertainties connected to the larger zone of avoidance. 

The initial conditions reconstructed from the NOG catalogue are found (analo- 
gously to those from the PSCz) to be consistent with Gaussian paradigm. We use the 
reconstructions to produce sets of initial conditions ready to be used for constrained 
simulations of our local Universe. 

Key words: cosmology: theory ~ galaxies: catalogues - large-scale structure of the 
Universe - galaxies: clustering - dark matter 



1 INTRODUCTION 

, In the gravitational instability scenario the observed large- 
' scale structure of the Universe forms by gravitational col- 
lapse of small perturbations imprinted in the early Universe 
by some mechanism like inflation. The perturbation field is 
believed to be a Gaussian stochastic process. In this case, 
all the cosmological information present in the perturbation 
field is given by its power spectrum, or by its Fourier trans- 
form, the two-point correlation function. On the other hand, 
gravitational evolution induces strong non-Gaussianities in 
the large-scale density distribution, easily recognizable both 
in the 1-point Probability Distribution Function (hereafter 
PDF) of the density, which is very similar to log-normal, and 
in the topology of large-scale structure, composed by voids, 
filaments and clusters (see, e.g., Sahni & Coles 1996). 

The specific pattern of primordial fiuctuations that has 
originated the structure observed in our Local Universe 
bears no cosmological information beyond that of the power 
spectrum, but a detailed knowledge of it is important for 
many purposes. Firstly, its reconstruction relies on the as- 
sumption of gravitational instability, but not on the Gaus- 
sianity of the initial field; the reconstructed density field can 



then in principle be used to test the Gaussianity of the ini- 
tial conditions (Nusser, Dekel & Yahil 1995, Monaco et al. 
2000). These test have not revealed any non-Gaussian fea- 
ture. However, they are not very sensitive to small signal, 
and the signal is degenerate with spurious non-Gaussian fea- 
tures induced by non-linear bias (Monaco et al. 2000). 

Secondly, the reconstruction can be used to run con- 
strained simulations that reproduce our local Universe. Such 
simulations are useful for at least two aims, namely to cal- 
ibrate distance indicators for non- homogeneous Malmquist 
bias (see, e.g., Strauss & Willick 1995; Kolatt et al., 1996), 
and to model galaxy formation in an environment that 
closely resembles our local Universe (Narayanan et al., 2001; 
Klypin et al. 2002; Mathis et al. 2002). Indeed, the local 
Universe is obviously the place where observations are most 
detailed, and where objects such as dwarf galaxies or LBS 
spirals are mostly observed. Any study of the environmental 
dependence of galaxy properties would greatly benefit from 
this kind of simulation. 

The density field at scales of from a few to a few tens 
of Mpc is evolving in the quasi-linear regime, in which 
the perturbations still retain memory of the initial condi- 
tions; however, the inversion of the gravitational evolution 
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in this regime is not trivial. This is due to many reasons: (i) 
the evolution has a significant degree of non-linearity; (ii) 
the presence of relaxed structure at small scales influences 
the reconstruction; (iii) the sampling of the density field by 
galaxy catalogues is not homogeneous, duo to the selection 
function, the presence of unobserved regions (like the Zone of 
Avoidance, hereafter ZOA) and to possible inhomogeneities 
of the catalogues; (iv) galaxies axe biased tracers of the to- 
tal density field, and the bias relation is poorly known; (v) 
distance estimates are unavailable for large complete galaxy 
catalogues, so that the density field is observed in the red- 
shift space. 

Standard N-body integrators are not suitable to solve 
the inversion problem, as any noise would be interpreted 
as decaying mode and amplified in the backward integra- 
tion. The proposed solvers are based either on the Zel'dovich 
(1970) approximation or on the minimum-action princi- 
ple (Pooblos 1989, 1990; Branchini & Carlbcrg, 1994). The 
most recent proposals arc PIZA (Croft & Gatzagnaga 1997), 
ZTRACE (Monaco & Efstathiou 1999), FAM (Nusser & 
Branchini 2000), MAK (Frisch et al. 2000); older proposals 
are reviewed in Narayanan & Croft (1999). These algorithms 
are able to invert the gravitational evolution on scales larger 
than 3-5 Mpc (Gaussian smoothing), where the highly- 
non linear regime starts to prevail. The inversion codes axe 
also used to find the density field in the real-space and the 
3D-field of peculiar velocities. 

The reconstruction has been attempted by many 
groups. Kolatt et al. (1996) used the IRAS 1.2 Jy redshift 
survey (Fisher et al., 1995), in order to recover the density 
field in real space, along with the velocity and the poten- 
tial fields. The initial conditions were then reconstructed 
by applying the BernouUi-Zel'dovich equation (Nusser & 
Dekel, 1992) to the peculiar velocity field. Then they forced 
the 1-point PDF of their resulting initial conditions to be 
Gaussian by applying a rank- and variance-conserving pro- 
cedure. Narayanan et al. (2001) applied the same approach 
to the analysis of the PSCz catalogue and tested their re- 
sults against simulations for various cosmologies and bias 
schemes. Mathis et al. (2002) started their analysis from 
the IRAS 1.2 Jy redshift survey. They smoothed heavily 
the density field and computed the velocity and the poten- 
tial fields. Through the application of Bernoulli-Zel'dovich 
equation they reconstructed the initial density field. They 
used this field as a constraint for the HofFmann-Ribak algo- 
rithm (Hoffmarm & Ribak, 1991), to obtain a realization of a 
Gaussian field that resembles the input one once smoothed 
on the same scale. Klypin et al. (2002) reconstructed the 
density and velocity fields applying the formalism of Wiener 
Filtering (see, Zaroubi et al, 1995) and the HofTmann-Ribak 
algorithm to the MARK III catalogue of peculiar velocities 
(Willick et al., 1997). Initial conditions were recovered by 
applying linear theory to the velocity field. 

In this paper wo present a new reconstruction of the lo- 
cal Universe which improves over the previous ones in many 
respect. We use the ZTRACE code (Monaco & Efstathiou 
1999; Monaco et al. 2000) to trace back in time the density 
field obtained from galaxy catalogues. This code is able to 
self-consistcntly invert the gravitational evolution preserv- 
ing the statistics of the initial density field, except in the 
high peaks which are dominated by highly non-linear dy- 
namics. The code has already been applied to the infra-red 



IRAS PSCz redshift catalogue (Saunders et al. 2000; Monaco 
et al. 2000), with the result of a Gaussian 1-point PDF of 
the reconstructed linear density field and no clear sign of 
distortion due to non-linear bias. 

For this reconstruction wo use both the PSCz and the 
optical NOG catalogue (Giuricin et al. 2000) to sample the 
galaxy density field. The advantage of this optical catalogue 
with respect to the infra-red one is an increase of density 
by a factor of 1.3 to 2 and a better relation of optical light 
with galaxy mass, which allows sampling of the early type, 
dust-poor galaxy population. The disadvantage is given by 
the larger ZOA and by the inhomogeneity of the sources, 
although the work of Marinoni et al. (1999) and Marinoni 
(2000) has demonstrated the very good degree of complete- 
ness of the sample. The use of two different catalogues allows 
us to test the robustness of the reconstruction to the sam- 
pling of the density field by different types of galaxies. 

The paper is organized as follows: section 2 describes 
the two catalogues used, and compares the galaxy density 
fields traced by them; the relative bias of the two catalogues 
is quantified, and this relation is used to fill the optical ZOA 
with the information given by the PSCz. Section 3 describes 
the tests performed on simulated galaxy catalogues to as- 
sess the accuracy of the reconstructed initial conditions in 
the Fourier space; in particular, we quantify the limit k- 
value at which the reconstruction is not reliable any more, 
and the corrections that are required to fix the flattening 
of the reconstructed peaks and to restore the power lost by 
smoothing. Section 4 shows the consistency of the initial con- 
ditions reconstructed from the NOG catalogue with Gaus- 
sianity, and the overall accuracy with which the structure of 
the local Universe is reconstructed by a set of constrained 
simulations. Finally, section 5 gives the conclusions. 

2 THE CATALOGUES 
2.1 PSCz and NOG 

Galaxies are biased tracers of the underlying density field, 
and this bias depends on the properties of the galaxy cat- 
alogue used. The PSCz catalogue (Saunders et al. 2000) 
meets the most important requirement, i.e. the homogene- 
ity of selection, in that it is extracted from the Point Source 
Catalogue of the IRAS satellite (Joint IRAS Science Work- 
ing Group 1988). In particular, all sources with 60/li fiuxes 
> 0.6 Jy which show galaxy-like FIR colors have been se- 
lected and re-observed in the optical. The final catalogue 
contains about 15000 confirmed galaxies over 84 per cent of 
the sky. Due to the homogeneity of the selection and the 
smallness of the ZOA (limited roughly to \h\ < 8° plus two 
thin stripes not observed by IRAS) the PSCz catalogue has 
been fruitfully used for many studies of large-scale structure 
(see Saunders et al., 2000 and references therein). Another 
advantage of the FIR selection relies in the high-luminosity 
cutoff of the selection function, which is weaker than the op- 
tical, and allows a sparse sampling to much larger distances. 
The main disadvantage relies in the poor relationship be- 
tween FIR light, mostly emitted by dust heated up by star 
formation or AGN, and galaxy mass. In other words, the 
FIR light is biased toward mid-type spirals (Marinoni, pri- 
vate communication), missing almost completely the early- 
type population. 
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Optical light is more tightly connected to galaxy mass, 
due to the presence of relations such as TuUy-Fischer or the 
fundamental plane. On the other hand, optical all-sky catar- 
logues must necessary rely on compilations of different sets 
of observations; this is a source of systematic error that may 
be difficult to subtract. Moreover, the ZOA in the optical 
is much wider than the IR. All-sky optical catalogues have 
been presented e. g. by TuUy (1988) and Santiago et al. 
(1995). We use the NOG catalogue (Giuricin et al. 2000), 
which has been selected from the LEDA database (see, e. 
g., Paturel et al., 1997) and then homogenized by Garcia et 
al. (1993), Marinoni et al. (1999) and Marinoni (2000). It is 
95 per cent complete to B = 14 and 98 per cent complete 
in redsliift. It is limited to |6| > 20° and cz < 6000 km s^^ , 
after whicli the selection function drops abruptly*. Giuricin 
et al. (2000) have extracted galaxy groups from the cata- 
logue using both a percolative, friends-of-friends algorithm 
(Huchra & GcUer, 1982) and a hierarchical one (Materne, 
1978). The catalogue has been used by Giuricin et al. (2001) 
for studies on redshift-space two-point correlation functions 
of galaxies and groups, by Giraxdi et al. (2002) for the anal- 
ysis of observational mass-to-light ratio of galaxy systems 
and by Marinoni, Hudson & Giuricin (2002) to obtain the 
optical luminosity function of virialized systems. 



2.2 The relative bias between PSCz and NOG 

It is important to quantify how the two catalogues map the 

galaxy density field, and in particular their relative bias. 
This is important to interpret possible differences in the re- 
construction of the initial conditions. Moreover, to obtain a 
full-sky coverage of the density field, it is important to fill 
in the optical ZOA with a galaxy field which resembles the 
true one as accurately as possible. The best way to do it is 
to use the information given by the PSCz itself, whose ZOA 
is much smaller. This requires an assessment of the relation 
between the galaxy density fields as described by the two 
catalogues. 

We can try to anticipate the result. The rms of galaxy 
counts in spheres of 8 Mpc, usually denoted as ag, 
has been measured to be ~ 0.9 for optical galaxies (Sza- 
lay et al., 2002) and ~ 0.8 for IRAS ones (Seaborne et al. 
1999, Sutherland et al. 1999). So, as long as both kinds 
of galaxies sample the density field with a roughly linear 
bias scheme, we expect to measure a relative bias parameter 

6rel — <78,opt/C8,FIR — 1.1. 

In the following we call $(L) the luminosity function 
of the catalogue. The presence of a flux limit fum deflnes 
a relation between the distance r and the absolute magni- 
tude M of a galaxy as the distance at which that galaxy is 
observed at the flux limit: 



imin(r) = 47rr^/lin 



(1) 



We define n(r) as the number density of objects in the fiux- 
limited catalogue at the distance r: 



n{> L, 



r+oo 
.in(r)) = / 



^{L)dL 



(2) 



* The cut in redshift is due in fact to the low degree of complete- 
ness of the LEDA database beyond that limit. 



Given a reference luminosity Ls, the selection function S(r) 
is defined as the fraction of galaxies more luminous than Lg 
that are visible at the distance r: 



S{t) = 



POO 

■J ma. 



x{Le,L™,„(r)) 



^{L)dL 



!ZHL)dL 



(3) 



while its inverse F{r) = 1/S{r) gives the number of galaxies 
more luminous than Lg present on average at distance r for 
each observed galaxy. For the NOG catalogue we use the 
luminosity function given by Marinoni (2000); in terms of a 
Schechter function: 



$(L) 



-L/L* 



(4) 



with parameters $* = 0.006 Mp(r'^{h/0.75f , M* = 
-20.6 + 51og(/i/0.75), a = -1.1. Ttie luminosity function 
in the infrared in given by Saunders et al. (1990). In this 
work we use the number density for the PSCz catalogue in 
the form (Saunders et al., 2000): 



n(> Lmin{r)) = 



(r/r* 



[1 + {r/r*)'<]Ph 



with parameters $* = 0.132 h^Mpc 



(5) 

* = 64.6 Mpc/h, 
a = —1.1, /3 = 3.90, 7 — 1.64. Figure 1 shows all these 
functions for the NOG and PSCz catalogues. The higher 
sampling rate of NOG is evident in this figure. 

The galaxy density fields of the two catalogues have 
been computed in different ways, so as to test the robust- 
ness of the result. The number of optical or IR galaxies has 
been computed in boxes which axe either sectors of spherical 
coronae or cubes. The first choice corresponds to a binning 
of galactic coordinates plus redshift; the binning of sky coor- 
dinates is shown in figure 2, while the binning in redshift is 
1200 km . The advantage of this choice is that all boxes 
are always either inside or outside the volume covered by 
both catalogues (i.e. that of NOG), and that it makes very 
easy to single out possible problems connected with inhomo- 
geneities in the selection, like differences between northern 
and southern hemispheres. The second choice corresponds 
to a binning of Cartesian coordinates, and allows to keep 
both the geometry and the volume of the boxes fixed. The 
catalogues have been immersed into a large box of side 120 
h^^ Mpc, with the Galaxy obviously at the centre. Then 
the large box has been subdivided into an integer number 
of small boxes. We consider only boxes that overlap for at 
least 75 per cent with the NOG volume. In the following we 
show the case of boxes of 10 ir^ Mpc. 
The density contrast is defined as: 



5 = 



A^a- 



(6) 



It has been computed with two different methods. In the first 
case, hereafter method I, we compare the number of galaxies 
TVgg^i found in the box at distance r with the expected one 
■^avei computed as: 



V 



/: 



$(L)dL 



(7) 



where $(i) is the luminosity function of the catalogue and 
imin(r) is defined in equation 1. In this case the shot-noise 
error associated with the density is: 
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Figure 2. The figure shows the binning in sky coordinates (galcLC- 
tic coordinates) employed in the analysis of relative bias between 
the catalogues. The thick line represent the celestial equator. Bin- 
ning in redshift coordinates is 1200 km s~^. 



e{l + S) 



gal 



(8) 



Alternatively (method II) we fix a limit luminosity Ls 
for the galaxies, and estimate the number of galaxies with 
L > Ls present in the box as: 



N. 



gal 



(9) 



The density contrast is obtained as in equation 6, where the 
average number of galaxies this time is: 



n!' = V 



'i>{L)dL 



(10) 



In this case the shot-noise error for the density is (Strauss 
& Willick 1995): 



e{l + 6) = 



Ni- 



1/2 



(11) 



The advantage of method I with respect to method II 
relics in the lack of necessity to define a limit luminosity (it 
is implicitly defined by the magnitude limit). In any case, 
we fix limit luminosities Ls (and correspondingly a distance 
Ts) for both catalogues for the following reason. In the NOG 
case, we know that dwarf galaxies trace a different density 
field with respect to the luminous ones, and that they are 
observed only in the inner part of the NOG volume, where 
they dominate in number. In this case, the density field in 
the inner part, which is oversampled, would be dominated 
by the dwarf population. This can induce some bias in the 
reconstruction. To avoid this wc fix a limit absolute magni- 
tude of -18.00 and cut both the catalogue and the luminosity 
function below this limit. Figure 3 shows the absolute mag- 
nitude - redshift plot for the NOG catalogue; the galaxies 
excluded are only 16.5 per cent, and all lie within Ts — 25.00 
h^^ Mpc. For the PSCz catalogue, the poor relationship 
between FIR luminosity and galaxy mass malces the previ- 
ous argument inapplicable. However, Rowan-Robinson et al. 
(1991) report an incompleteness of the PSCz catalogue for 
the nearest galaxies; following their suggestion we set Ts = 6 
h'^ Mpc. 

Figures 4 and 5 show the resulting correlation between 
PSCz and NOG densities, for the two kinds of boxes and 
the two methods I and II. Each panel shows the bisector line 
and the standard rank correlation coefficient r for the cor- 
relation. We show the correlation both in 5 and in log(l-|-(5) 
because it is more convenient to describe a bias relation in 
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Figure 3. Absolute Magnitude - Redshift plot for the NOG cat- 
alogue. Objects with absolute magnitude greater than —18 are 
excluded from the analysis. They represent 15 % of the total cat- 
alogue. 



Figure 5. Correlation between PSCz and NOG density fields 
in a bilogaritmic plane. Methods and volumes employed in the 
analysis are indicated at the top of each panel. Rank correlation 
coefficients are reported at the bottom of each panel. 
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Figure 4. Correlation between PSCz and NOG density fields in 
a linear plane. Methods and volumes employed in the analysis are 
indicated at the top of each panel. Rank correlation coefficients 
are reported at the bottom of each panel. 



the log(l + S) variable, as the 5 > — 1 constraint is auto- 
matically satisfied. For all four methods the correlation is 
very good, with correlation coefficients exceeding 0.65 in all 
cases, and as high as 0.8-0.9 in the case of Cartesian boxes. 
The greater scatter found for spherical volumes with respect 
to the cartesian ones is due to their different shapes (see 
section .3.f, below). This effect could be induced in prin- 
ciple by a scale-dependent bias. However, we do not find 
significant variations in the bias parameter when estimated 
within cartesian boxes of size from 7.5 to 2()h — IMpc. Wo 
note that this scale range covers most of the range probed 
by the spherical volumes. The quality of fitting is stable 
with distance and position on the sky, suggesting that the 
inhomogoncity of NOG induces no significant error in the 
reconstruction of the density (sec also Marinoni 2000). 

All the plots have been fit by a linear relation, to deter- 
mine the bias parameters. In particular: 



Onog 

log(l-|- 5nog) 



= aiog -I- 6iog log(l-l- 5ps, 



(12) 



It is worth stressing that the two aun and 6iin or aiog and 
6iog parameters are not independent, as they are subject to 
the constraint (S) = for both density fields. Moreover, the 
linear coefficients can be obtained by the first-order term of 
the Taylor expansion of the log relation. Table 1 gives the 
a and b coefficients for all the regressions. The b bias coeffi- 
cients are all consistent, and range from 1.10 to 1.20, with an 
error of 0.1, with only one exception, corresponding to a low 
r-value for the correlation. The no-bias CEise is never signifi- 
cantly rejected. We conclude that the bias relation between 
the FIR and the optical galaxy density fields is: 



brel = 1.1 ± 0.1 



(13) 
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Table 1. a and b coefficients for all the regressions (equations 12). 
The fits are performed assuming Poisson errors on both axes. 
With the simbol ar we indicate the scatter of 5nog around the fit. 



Binning 5 hii„ aii„ Cr 

Sky coord. Met. I 1.08 ± 0.10 0.04 ± 0.06 1.31 

Sky coord. Met. II 1.20 ± 0.10 0.17 ± 0.07 1.27 

Cartesian coord. Met. I 1.17 ± 0.10 0.10 ± 0.06 0.35 

Cartesian coord. Met. II 1.19 ± 0.10 0.13 ± 0.07 0.37 



Binning <5 hiog aiog Cr 

Sky coord. Met. I 0.94 ± 0.09 -0.01 ± 0.05 0.36 

Sky coord. Met. II 1.11 ± 0.11 0.10 ± 0.06 0.38 

Cartesian coord. Met. I 1.16 ± 0.10 0.02 ± 0.06 0.23 

Cartesian coord. Met. II 1.12 ± 0.10 0.06 ± 0.06 0.23 



This is consistent with our expectation of broi — 1.1. 

We notice that the points that lie most significantly out 
of the correlation are those contained in the first bin in red- 
shift (cz < 1200 km If we exclude these points from 
the analysis we find no relevant variation in the values of fit 
parameters, while the rank correlation coefficient increases 
significantly. It is possible that this behaviour is related to 
the incompleteness of PSCz in the nearest region or by the 
presence of the E-rich Virgo cluster region. We have also 
tried to change the number and shape of the volumes con- 
sidered, obtaining a set of values for the a and & coefficients 
consistent with those given in table 1. 

An estimate of the statistics for these regressions 
gives in all cases a high value, with a very low associ- 
ated probability. We interpret this as a signature of relative 
stochastic bias: the two density fields have an intrinsic scat- 
ter beyond that induced by sparse sampling. We arbitrarily 
try to double the Poisson error given by equations 8 and 11, 
and obtain acceptable values. This means that the in- 
trinsic scatter of the two density fields is similar to the shot 
noise of the reconstruction. In order to quantify the degree 
of stochastic bias between the distribution of optical and in- 
frared galaxies, we use the coefficient Oh = ar/os, where ar 
is the scatter of 5nog around the best fit line (see Table 1) 
and as is the variance of the (5pscz- We obtain: 

a6~0.25. (14) 

The same parameter was used by Dekel & Lahav (1999) to 
quantify the stochastic bias of galaxies with respect to dark 

matter. Interestingly, our value is similar to the prediction of 
Somerville ct al. (2001) for the stochasticity of optical galax- 
ies with respect to dark matter, which results in their case in 
the range 0.18 — 0.35 (depending on the absolute magnitude 
cut), with a stronger signal from optically bright galaxies. 
Of course this result strongly relies on the hypothesis that 
Poisson errors are correct. Furthemore, non-linearity in the 
bias relation (that is not observed in our data) would lead 
to an overestimate of Ob- 



2.3 Filling the optical ZOA and the external 
regions 

Due to the non-local character of gravity, it is necessary to 
fill in the ZOA with a mock galaxy distribution. Indeed, if 
the ZOA were naively treated as a void, it would push all 
the surrounding matter aside, thus biasing heavily the re- 
construction. Filling it with random galaxies, thus neglect- 
ing completely the structure present, would also distort the 
reconstruction. To fill in the FIR ZOA we have used the 
scheme of Branchini et al. (1999), also used in Monaco et 
al. (2000), which consists in interpolating the density field 
of two stripes above and below the ZOA. 

In the case of NOG, the optical ZOA is so wide that 
such an interpolation would result in a poor representation 
of the true density field. Moreover, NOG is limited to cz < 
6000 km s^^, while PSCz sparsely samples the density field 
to a distance of ^ 200 — 300 Mpc, a region which gives 
an important contribution to local bulk motion and tides. 

It is convenient to use the information provided by the 
PSCz catalogue, together with the observed relation of op- 
tical and FIR galaxy density fields, to generate a mock cat- 
alogue of optical galaxies that fills the ZOA and extends to 
distances larger than 60 Mpc. This is done as follows. 
The PSCz catalogue has been immersed into a box of side 
250 Mpc. Its density field has been computed on a grid 
of 128^ points by using the adaptive smoothing scheme al- 
ready used by Canavezes et al. (1998), Monaco & Efstathiou 
(1999) and Monaco et al. (2000). In this scheme a Gaussian 
smoothing kernel is associated to each galaxy, with a refer- 
ence radius i?rcf that satisfies the relation: 

— il?ef(r)n(> Ln,in(r)) = 10 (15) 

(so as to have on average 10 galaxies within a distance Rid). 
Next, an adaptively refined radius i?adap is computed, such 
as to contain exactly 10 galaxies within itself. Then, the 

density field is computed by summing over all kernels. In 
this way the density field is computed with constant signal- 
to-noisc. 

After this, the volume of the box is divided into small 
cells, which are populated by galaxies according to the prob- 
ability: 

P = V ■ Snog ■ (1 + 5nog) f^^.-. 
= V ■ Snog ■ [(1 + (^pscz)"-' 10"">] ^ 

The boxes are kept small enough so that the probability of 
hosting a galaxy is always smaller than one: a volume of 
about 1 square degree in sky coordinates and 50 km in 
redshift space is small enough. Once the box is randomly 
selected to contain a galaxy, the object is put into a ran- 
dom position within the box. With this procedure the mock 
catalogue follows both the NOG selection function and the 
PSCz density field, rescaled to the NOG one through the 
relations defined above. 

In practice, we do not use the a and b parameters found 
by the linear regressions. The 6rei parameter is set either to 
1 or to 1.2, to be able to test the robustness of the results 
to uncertainties on the relative bias. The corresponding Orel 
parameters are chosen by requiring a good fit of the selec- 
tion function for ten different realizations of the filling mock 
catalogue. In this way we obtain Orei = —0.03 for 6rei = 1-2, 
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Figure 6. Completions of the NOG catalogue. Left panel shows the NOG catalogue with the ZOA filled with the information from PSCz 

catalogue in the hypothesis that the relative bias parameter 6[.cl is 1- Right panel shows the NOG catalogue with the ZOA filled with 
the information from PSCz catalogue in the hypothesis that the relative bias parameter by-gi is 1.2. 



while Orel = 1 for 6rei = 0. We do not account for stochas- 
ticity in the filUng procedure. 

Figures 6 shows the original NOG catalogue completed 
with two mock catalogues extracted from the PSCz one, 
assuming 6rei = 1 or 1.2. The structures visible are continued 
into the ZOA, although the voids are not as neat as in the 
real catalogue. Moreover, the structures in the 6rei = 1-2 
case show a slightly higher contrast. 



3 TESTS ON SIMULATED CATALOGUES 

The ability of ZTRACE to reconstruct the initial conditions 
of our local Universe was demonstrated by Monaco & Ef- 
stathiou (2000). In the following we perform further tests 
aimed to quantify the accuracy gained by using the denser 
NOG catalogue with the ZOA filled by the PSCz and to cal- 
ibrate the corrections for peak flattening and smoothing. In 
particular, we test the goodness of the reconstruction of the 
initial density field in the Fourier space, considering not only 
the moduli of the modes (whose average is the power spec- 
trum) but also their phases, which are mostly responsible 
for shaping the observed large-scale structure. 

3.1 Simulated catalogues 

For these tests we perform the ZTRACE reconstruction on 
galaxy catalogues extracted from a cosmological N-body 
simulation. The simulation has been run using the Tree grav- 
itational part of the GADGET code^ (Springel, Yoshida 
& White 2001). We simulate a ACDM model with Qrn = 
1 - Qa = 0.3, Qbar = 0.019 ft"^, h = 0.7 and as = 0.8, 
using 256^ dark matter particles within a box having co- 
moving size of 250h~^Mpc. The Plummcr-cquivalent soft- 
ening length for the computation of the gravitational force 
has been set to ep; = 30/i~^kpc fixed in comoving units. 
Initial conditions have been generated at Zi ~ 16 using the 
GRAPHIC2 package by Bcrtschinger (2001). 

Mock catalogues are extracted following the selection 
function either of the NOG or the PSCz catalogue, and a 

t http: //www.MPA-Garching.MPG.DE/gadget/ 



bias scheme analogous to that described in section 2 (see 
equation 16) , where the probability that a particle is selected 

as a mock galaxy is: 

Tlmax 

Here rimax is the mean density of particles in the simulation 
and: 

P6 = a.in,(l-h<5)'''='--'^ (18) 

In this case the density field 6 is the final one traced by 

all the particles in the simulation, cic-interpolated on a grid 
of 128"' and smoothed over 3 h^^ Mpc. We have checked 
that with our cosmology the variance of particle counts on 
boxes of side L is equal to the variance of the density field 
Gaussian smoothed on scale R when L ~ 3.57?; so 3 Mpc 
Gaussian smoothing is roughly equivalent to the 10 h^^ Mpc 
boxes shown in figure 4 and 5. With this bias scheme mock 
galaxies are power-law biased with respect to the dark mat- 
ter. Again, the 6sim parameter is fixed to some value, while 
the asim parameter is chosen so as to reproduce the selection 
function. 

The volume of the NOG catalogue (a sphere of 60 
Mpc) is only 6 per cent of the volume of the simulation, 
so that it is safe to extract ten independent catalogues from 
the simulation. We extract four sets of ten catalogues, all 
centred on the same set of ten positions (periodic boundary 
conditions axe assumed). Two of them are extracted with the 
NOG selection function, one unbiased (6sim = 1, ftsim = 0) 
and the other with bsim ~ 1.2, Osim = —0.03; two more are 
selected with the PSCz selection function and either unbi- 
ased or with 6sim = 0.8 and fisim = 0.03. 

After defining a galactic plane (the x-y plane defined by 
the box) the optical ZOA and the part with d > 60 Mpc 
arc taken out from the NOG catalogues, and then filled with 
the corresponding simulated PSCz catalogues following the 
procedure already outlined in section 2.3: (i) wc compute 
the density of the corresponding simulated PSCz catalogues 
by adaptivcly smoothing with fixed signal-to-noise; (ii) we 
extract a mock catalogue which samples the PSCz density 
with a given bias scheme; (iii) with it we fill the optical 
ZOA and the external parts. The real PSCz is also subject 
to a procedure to fill the ZOA, but we neglect to apply this 
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Table 2. Simulated Catalogues 



Catalogue 


^sim 


NOG unbiased 


1 


NOG biased 


1.2 


PSCz unbiased 


1 


PSCz biased 


0.8 


Catalogue 





NOG unbiased filled with PSCz unbiased 1 
NOG biased filled with PSCz unbiased 1.2 
NOG unbiased filled with PSCz biased 1.2 



to the volume geometry. The outcoming slopes of the re- 
lations are for the four cases brci ~ 1.22 ± 0.06 (Method 
I, Spherical), b^ei = 1-18 ± 0.05 (Method II, Spherical), 
6rei = 1.11 ±0.06 (Method I, Cartesian) and 6rei = 1.06±0.03 
(Method II, Cartesian). This leads to an average estimate 
of 6rei = 1-14 ± 0.06, consistent with the true 6rei = 1-2 
value, although we notice (also in the other realizations) a 
slight tendency to underestimation of 6rei- Moreover the 
tests give reasonable values with high associated probabil- 
ity, confirming that the scatter in this case is due to Poisson 
errors; this supports our interpretation of the scatter in the 
data as stochasticity of bias, although with all the caveats 
highlighted above. 




-2-1012 -2-1012 



Log(l+(5 ) Log(l+(5 ) 

Figure 7. As in figure 5 for one simulated pair of catalogues; 
NOG with bsim = 1-2, PSCz unbiased. 

to the mock catalogues. Indeed, it was shown in Monaco 
& Efstathiou (1999) that the impact of this filling proce- 
dure is modest, and surely smaller than the corresponding 
procedure for the NOG. The relative bias schemes assumed 
for the filling procedure are an unbiased one, applied to the 
pairs of unbiased NOG and PSCz catalogues, and a biased 
one, applied either to the biased NOG and unbiased PSCz, 
or to the unbiased NOG and (anti)biased PSCz. In total, we 
reconstruct seven sets of ten simulated catalogues; these are 
listed in table 2. 

As a first check wo repeat the procedure described in 
section 2.2, i.e. we divide the volume into subvolumes and 
compute the "galaxy" density contrasts for pairs of simu- 
lated catalogues. Figure 7 shows the bilogarithmic corre- 
lation of density contrasts for one realization of unbiased 
PSCz and biased (bsim = 1-2) NOG. As in figure 5, spherical 
volumes show a greater scatter than cartesian ones, espe- 
cially when Method I is used. As anticipated above, this 
confirms that the variation in scatter is mostly connected 



3.2 ZTRACE reconstruction 

Here we outline briefly the main steps of the ZTRACE re- 
construction. This applies to the simulated catalogues as 
well as to the real ones. Full details are given in Monaco & 
Efstathiou (1999). 

First the relaxed groups are collapsed, so as to remove 
the "finger of God" features present in the redshift space; 
this is known to improve the reconstruction (Gramann et al. 
1994). In all cases but the real NOG catalogue, to select the 
groups we use a standard friends-of-friends code with linking 
lengths 0.5 and 3 Mpc respectively tangential and along 
the line of sight. In the NOG case we use the groups found 
by Giuricin et al. (2000). This operation is obviously carried 
over the catalogues, and not over the filling of the ZOA. 

Second, the catalogues are smoothed with the same 
adaptive scheme as that described in section 2.3. This time 
the reference smoothing radius is kept constant to -Rmax 
within a distance dmax, such as to have 47r/3 -Rmax'^(> 
imin(dmax)) = 10, and theu changed as in equation 15 so 
as to keep the signal-to- noise constant. The adaptive refine- 
ment is done so as to have 47r/3 R^i,fn{> Lmin{r)) galax- 
ies within -Radap- In this way the inner part (r < dmax) is 
adaptively smoothed with constant reference smoothing ra- 
dius, while the outer parts are smoothed with a reference 
smoothing radius which increases with distance. For choices 
of dmax = 30 or 60 Mpc we have -Rmax = 3.17 and 5.57 
h-^ Mpc for NOG, 3.83 and 5.99 h'^ Mpc for PSCz. 

Third, the density field is given as an input to the 
ZTRACE code, thus obtaining a reconstruction of the initial 
conditions. These are given in terms of the quantity: 



where (5(x, U) is the density contrast field at a very early time 
ti, when linear theory holds, and D(ti) is the linear growing 
factor at ti. This quantity, called linear density contrast, is 
constant in linear theory, and is equal to the linear extrap- 
olation of the initial conditions to z = Q (where D(t) = 1 
by normalization). At variance with Monaco & Efstathiou 
(1999), the initial conditions are reconstructed on a grid of 
a 128^ (instead oi M^)^ . 

■t This way the smoothing at the centre of coordinates, necessary 
to ensure convergence of the iteration, affects only the density 
field within ~ 5 Mpc. 
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Table 3. Parameters of the correction for the flattening of the 
pealis 



Catalogue 


dmax 


= 60Mpc/h 


C^max 


= 30Mpc/h 




CI 


C2 


CI 


C2 


NOG unbiased 


1.03 


6.50 


0.95 


4.20 


NOG biased 


1.00 


4.70 


0.90 


3.70 


PSCz unbiased 


1.05 


4.80 


0.97 


4.95 


PSCz biased 


1.05 


6.00 


1.00 


4.70 


NOG unbiased filled 


1.10 


5.10 


1.10 


3.50 


with PSCz unbiased 










NOG unbiased filled 


1.04 


5.30 


0.95 


3.90 


with PSCz biased 










NOG biased filled 


1.0.5 


4.00 


0.95 


4.00 


with i*S('z uiil)ias(;(l 











3.3 Fixing the flattening of the peaks 

The ZTRACE algorithm is known to preserve the statis- 
tics of the initial conditions for density contrasts lower than 
5; ~ 1, but it flattens the high-density peaks. This is due 

to the fact that high peaks arc sites of multi-streaming, a 
regime in which any reconstruction based on the mildly non- 
linear evolution fails; ZTRACE assumes that the density 
field is always in single-stream regimes. This distortion in- 
duces some degree of non-Gaussianity in the reconstruction 
that must be removed. 

To achieve this we compute the PDF of the recon- 
structed initial conditions for all the ten reconstructions of 
a given set of simulated catalogues, and from these the aver- 
age initial PDF. Galaxy bias is known to induce distortions 
in the initial PDF (Monaco et al. 2000), but the effect is neg- 
ligible for the moderate bias schemes used here. To decrease 
the scatter in the PDF from one realization to the other, 
we rescale all the PDFs to zero mean and unit variance, use 
the median instead of the mean and compute the variance 
using only the underdensities, were the density contrast is 
best reproduced. The reconstructed initial PDF is found to 
be always consistent with Gaussian for 5 < 0. For 6 > we 
choose a parametric transformation 5' = f{S) such that the 
PDF of the 6' variable is Gaussian to a good approxima- 
tion. For this transformation we use a branch of hyperbola, 
subject to the constraint of vanishing at 5 = with first 
derivative equal to unity (so as to join smoothly the bisec- 
tor line), and to have a flat asymptote at (5 — > oo: 

J ax — \/ c?a;2 -|- C2(ci — l)a; -|- Co -\- ci if a; > 

iix<0 (20) 

The ci and C2 parameters determine respectively the cur- 
vature of the function and the level of the asymptote. The 
parameters have been found for each reconstructed PDF; 
they are listed in table 3. 

Figure 8 shows, as an example, the average recon- 
structed PDF for the unbiased NOG catalogues with the 
ZOA fllled with the unbiased PSCz density flelds, before 
and after the fix for peak flattening. Figure 9 shows the 
scatter-plot of true and reconstructed density for one of the 
simulated catalogues, with the best-fit transformation su- 
perposed. 




Asim 

Figure 9. Scattcr-plot of true and reconstructed linear densities 
for one of the simulated catalogues (NOG unbiased filled with 
PSCz unbiased), with the best-fit transformation. 

Finally, we notice that our procedure is conceptually 
different from the Gaussianization technique of Weinberg 
(1992), where a Gaussianity of the PDF is directly forced on 
a linearly reconstructed density field. We apply our "Gaus- 
sianization" procedure only on the high-density peaks (the 
5 < density field is completely untouched) to remove the 
already analized effect of peak flattening. 

3.4 Accuracy of the reconstruction in the Fourier 
space 

Due to the inliomogeneous nature of the smoothing per- 
formed on the original catalogues, the initial conditions are 
not recovered down to a scale which is constant in space. 
Then, to test the reconstruction in the Fourier space with- 
out mixing information from different scales we pad to zero 
both the true and the reconstructed linear density flelds be- 
yond dmax- We also pad in some cases the ZOA, to test the 
reconstruction in the region which is not subject to the fill- 
ing procedure. This padding of course induces some phase 
correlation at large scales. To be sure that any positive cor- 
relation of moduli and phases that we find is due to a true 
signal, we compare the reconstructed linear density field also 
with the true one centred in a point which is half the box size 
away from the centre of the simulated catalogue considered. 
In this case the reconstructed and true initial conditions 
are completely independent, and any correlation must be 
ascribed to the padding procedure. These initial conditions 
will be referred to as "random" in the following. 

The true, reconstructed and random initial conditions 
have been FFT-transformed, and the modes Sk have been 
compared one by one in terms of modulus |<5k| and phase 
ifk. To quantify the agreement the modes have been binned 
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Figure 8. Average reconstructed PDFs. Left panel shows the average reconstructed PDF for "NOG- unbiased filled with PSCz-unbiased" 
simulated catalogues. The non-Gaussianity is due to the flattening of high-density peaks. Right panel shows the same PDFs after the 
correction of peak flattening (with coefficients ci = 1.10, C2 = 5.10; see table 1). In both panels dotted lines show a Gaussian of zero 
mean and unit variance. 



into intervals of k, and for each bin the following statistics 
have been computed: 



odule 



|5j.cc|2 



(21) 



(22) 



Figure 10 shows as an example these statistics for the 
unbiased NOG catalogues filled with unbiased PSCz's. The 
curves arc averages of the statistics over the ten realizations, 
the 1-0" scatter is also shown. Cross points refer to the com- 
parison of true and reconstructed fields, square points to the 
comparison of random and reconstructed ones. The thick 
lines show the expectation for two completely uncorrelated 
fields, which for ^module is: 



"module 



psim(J,)2 

^l-exp"*^^^j +exp~ 

(the last passage is done assuming P'^"" ~ ps'm gxp 
which is approximately true) . For the phases we obtain: 

rfuncor 
^phase 

~ 103.92 



(23) 



-k^R-' 



s^i:'^: = jd4>ij #2 5^^(01 - 02)' 



(24) 



The integral is computed taking into account that — Tt < 

(01 — 02) < TT. 

The first statistics quantifies the correlation of the fiuc- 
tuations of the modules of the modes around the average 



value (i.e. the power spectrum), the second statistics quan- 
tifies the correlation of phases. From figure 10 it is clear that 
both phases and fluctuations of modules are significantly re- 
constructed at large scales, while the correlation tends to the 
random value at large fe-values. The transition happens at 
a fc- value similar to 1/R, so that in this case the worsen- 
ing of the agreement is mostly due to smoothing. Finally, 
phases are apparently reconstructed more accurately than 
the fluctuations of the modules. 

It is important to quantify the fc- value at which the cor- 
relation is lost. We have tried many possible definitions for 
this fciim quantity, mostly based on the agreement between 
the phases, which we consider more important than the fluc- 
tuations of the modules. Namely, we have defined kum as the 
value at which the correlated curve is for the first time within 
(or for the last time beyond) 2a of the uncorrelated curve, 
with a that relative of either of the two curves, or the value 
at which the average correlated curve is at a fixed fraction 
(say 90 per cent) of the value relative to uncorrelated fields 
(equation 24). All these definitions give similar values for 
fciim. In particular, the last definition is the most stable one, 
as it is not based on the noisy curves relative to the scatter. 
In the following only the results based on the last definition 
will be presented. 

Table 4 and 5 give the resulting fciim for the seven sets 
of simulated catalogues, the two choices of c?max and i?max 
and the two cases of no padding or padding of the ZOA. 
For reference, we give also the quantity feum^^max, which is 
a good indicator of the goodness of the reconstruction. We 
conclude that: (i) for dmax = 60 Mpc the reconstruc- 
tion is valid up to feiim^imax — 2, i.e. at scales as small as 
half of the (Gaussian) smoothing radius; (ii) in the ideal 
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Figure 10. Statistic of reconstruction. The left jjaiicl siiows Smodule- The right panel shows Sphase- The curves are averages of the 
statistics over the ten reahzations. Cross points refer to the comparison of true and reconstructed fields, square points to the comparison 
of random and reconstructed ones. Continuous lines represent the l-cr scatter from the average value for the comparison of true and 
reconstructed fields, dashed lines for the comparison of random and reconstructed ones. The thick lines show the expectation value of 
the statistics. 



Table 4. Results for kum; dmax = dOMpc/h 



Catalogues 






^'lim-^niax 


ZOA i)a<l<lc<l 


^'lim-^max 

ZOA i)a<l<lc<l 


PSCz unbiased 


5.99 


0.36 


2.16 






PSCz biased 


5.99 


0.37 


2.22 






NOG unbiased 


5.57 


0.38 


2.12 






NOG biased 


5.57 


0.38 


2.12 






NOG unbiased filled 
with PSCz unbiased 


5.57 


0.37 


2.06 


0.37 


2.06 


NOG unbiased filled 
with PSCz biased 


5.57 


0.36 


2.01 


0.37 


2.06 


NOG biased filled 
with PSCz unbiased 


5.57 


0.34 


1.89 


0.36 


2.01 



case of no ZOA and no bias, the denser sampling of NOG 
with respect to PSCz gives a tiny improvement to the re- 
construction in terms of feum, with a decreasing fcumiimax; 
(iii) the correction of the ZOA makes the advantage of NOG 
almost vanish; (iv) the gain in fenm obtained by smoothing 
at a smaller scale (with dmax — 30 Mpc) is modest, cor- 
responding to a decrease of fcum^^max, which confirms that 
scales smaller than ~ 5 /i"'^ Mpc are dominated by highly 
non-linear dynamics; (v) in the case of biased catalogues the 
reconstruction works to a smaller value, because the on- 



set of non-linearity is anticipated, while the opposite is true 
in the case of antibiased catalogues^ . 



S It is fair to recall here that ZTRACE works always in the hy- 
pothesis of no bias, so that a biased catalogue gives a more non- 
linear density field. 
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Table 5. Results for kum', dmax = 30Mpc/h 



Catalogues 


RmB^(Mpc/h 


hi^{h/Mpc) 


^lim-^max 


himih/Mpc) 
ZOA padded 


ZOA padded 


PSCz unbiased 


3.83 


0.38 


1.45 


- 


- 


-T ov^z Diasca 


o.oo 


n /in 

U.4U 


l.Oo 






NOG unbiased 


3.17 


0.41 


1..30 


- 


- 


NOG biased 


3.17 


0.41 


1.30 






NOG unbiased filled 
with PSCz unbiased 


3.17 


0.38 


1.20 


0.40 


1.27 


NOG unbiased filled 
with PSCz biased 


3.17 


0.38 


1.20 


0.38 


1.20 


NOG biased filled 
with PSCz unbiased 


3.17 


0.38 


1.20 


0.38 


1.20 



by the reconstruction. The same conclusions hold for the 
other cases. 



0.2 
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Figure 11. Distribution of the phases of true (dotted lines) and 
reconstructed (continuous lines) initial conditions for k < fciim- 
The distributions are averages over 10 simulated catalogues. We 
show the case for unbiased NOG catalogues filled with unbiased 
PSCz. 



3.5 The distribution of the phases 

It is useful to check the Gaussianity of the reconstructed ini- 
tial conditions directly in terms of distribution of the phases. 

Figure 11 shows, for the unbiased NOG catalogues filled 
with unbiased PSCz, the average distribution of phases for 
k < fciim, together with the average distribution of the 
phases of the true initial conditions padded beyond dmax- 
We have checked that both distributions arc consistent with 
flat by FFT-transforming them and comparing their spec- 
trum with a white-noise one. We conclude that the flatness 
of the distribution of the reconstructed phases is preserved 



3.6 Correction of the povirer spectrum for 
smoothing 

Figure 12 shows the power spectra of the padded true, 
padded true smoothed (over 7?max) and reconstructed initial 
density fields averaged over the ten realizations in the case of 
unbiased NOG catalogues filled with unbiased PSCz. The re- 
constructed power spectrum is similar to the smoothed one, 
confirming that the adaptive smoothing scheme adopted by 
ZTRACE is very similar to smoothing in the (Lagrangian) 
space of initial conditions. Similar results are obtained in 
other cases. However, the small discrepancies visible in the 
figure and present also in the other cases are sometimes 
larger than the typical error associated to the power spec- 
trum (i.e. the rms over the ten realizations). 

To restore the power lost by smoothing, it is possi- 
ble as a first approximation to simply multiply the FFT- 
transformed field by the inverse of the Gaussian smoothing 
kernel, i.e. by exp[(fcJ?max)'^/2]. This turns out to be a good 
approximation for the two NOG catalogues. In the PSCz 
case the correction is improved by multiplying the Gaussian 
kernel by the ratio of the reconstructed and true (smoothed) 
padded power spectra, as those shown in figure 12. To check 
the accuracy of these corrections we discard from the FFT- 
transformed reconstructions all the modes with k > feum- 
The density fields are then completed with Gaussian ran- 
dom fields with power spectra truncated below fcum- The 
so-completed fields are then transformed them back to the 
real space, padded beyond dmax and FFT-transformed again 
to compute the power spectra. We obtain smooth spectra in 
all the three cases, without any jump at kum- 



4 RESULTS 

The three catalogues (PSCz and the two NOG comple- 
tions with either ferei = 1 or 1.2) have been processed as 
described in section 3; for sake of clarity we sum up the 
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Figure 12. Power spectra. Dashed line refers to the spectrum of 
true padded initial density field; dotted line refers to the spec- 
trum of true padded initial density field smoothed over 6 Mpc/h; 
continuous line refers to the average spectrum of reconstructed 
initial density fields over 10 simulated catalogues. We show the 
case for unbiased NOG catalogues filled with unbiased PSCz. 

procedure again, (i) The relaxed groups have been col- 
lapsed to spheroids; (ii) the catalogues have been adaptively 
smoothed, with reference smoothing radius kept constant 
within either 60 or 30 Mpc; (iii) the density fields have 
been processed by ZTRACE to generate the initial condi- 
tions; (iv) the flattening of the peaks has been fixed as de- 
scribed in section 3.3; (v) the power lost by smoothing has 
been restored as described in section 3.6. 

4.1 The PDF of the initial conditions from the 
NOG catalogue 

It is interesting to check the Gaussianity of the initial con- 
ditions reconstructed from the NOG catalogue, as done by 
Monaco et al. (2000) for the PSCz. 

Figure 13 shows the reconstructed PDFs for the two 
NOG catalogues (filled assuming either 6rei = 1 or 1.2) for 
rfmax ~ 30 or 60 Mpc; at this stage the flattening of 
the peaks and the restoration of power lost by smoothing 
have not been applied yet. The curves are compared with 
the average PDFs computed from the simulated catalogues. 
Consistently with Monaco et al. (2000) , no distortion of the 
PDF is detected beyond that induced by ZTRACE. There is 
a small discrepancy at 5 ~ 0, but it is valid only for one or at 
most two points, and does not resemble any of the patterns 
induced by realistic bias schemes (see Monaco et al. 2000); 
we do not consider this feature as significant. 

Finally, figure 14 shows the phase distribution of the two 
sets of initial conditions padded to zero beyond dmax = 60 

Mpc, for k < feum, where fcum is taJfen from the corre- 
sponding case of table 4. For sake of comparison we show 



also the phase distribution of the initial conditions padded 
beyond dmax for one of the simulated catalogue centres. It 
is apparent that the weak difference from flatness of the re- 
constructed phases is consistent with it being an effect of 
padding (see also figure 11). A similar result holds for the 
PSCz reconstruction. 

4.2 Simulations of the local Universe 

The three reconstructions of the initial conditions of our lo- 
cal Universe have been fixed for flattened peaks, as described 
in section 3.3, and for smoothing, as described in section 3.6. 
They have then been FFT-transformcd and truncated at 
feiim. Finally, the modes at fc > feum have been extracted as- 
suming a ACDM power spectrum with the parameters given 
in section 3.1 and random phases. 

The so obtained density field can be used to generate 
displacements for a set of particles, so as to produce ini- 
tial conditions for a constrained N-body simulation of our 
local Universe. A more refined way to do it relies on the 
Hoffmann & Ribak (1991, 1992) approach to generate con- 
strained realizations. This would produce a Gaussian density 
field that is not oversmoothed at large distances. However, 
in a multi-grid simulation the density field at large distance 
would anyway be sampled by rather massive particles, thus 
removing the intermediate scale recovered by the HofFmann- 
Riback method. 

A comparison of a full-blown simulation with the prop- 
erties of our local Universe, like the analysis of Narayanan 
et al. (2001), would require a careful assessment of the bias 
relation between dark matter and galaxies, and is consid- 
ered beyond the scope of this paper. However, to perform a 
"quick and dirty" check of the goodness of the reconstruc- 
tion, we follow the evolution of 30 different realizations of 
the initial field (ten for each different reconstructed cata^ 
logue) using the PINOCCHIO tool, recently proposed by 
Monaco et al. (2002)^. This gives a very fast and accurate 
approximation of the final distribution of dark matter halos 
in terms of masses, positions, velocities and merger histo- 
ries, and has the advantage of running in a very short time 
compared to a standard N-body code. We are interested 
in the reconstruction of the large-scale structure, so high 
resolution is not required for this test. For this reason we 
run PINOCCHIO on 128'^ grids, for a mass resolution of 
~ 10^^ Mo (/i = 0.7). In this case, the smallest reliable halo 
reconstructed by PINOCCHIO (as well as by an equivalent 
simulation) is as massive as ~ 3 — 5 x 10"'^'' Mq\ this mass 
scale is still dominated by the small-scale power which is 
randomly added to the ZTRACE reconstruction. 

Figure 15 shows the reconstruction of the NOG volume 
along the super-galactic plane. Upper panels show 3 real- 
izations based on the NOG catalogue with brd = 1.2; mid 
panels show the same for brei = 1.2; lower panels show the 
same for the PSCz catalogue. Stars indicate the positions of 
halos reconstructed by PINOCCHIO, with mass grater than 
8 X 10^^. For sake of comparison we report in figure 16 the 
location of NOG galaxies. The main structures observed in 
the supergalactic plane are correctly reconstructed. In par- 
ticular, we find that the locations of the main neighbour 

^ http://www.daut.univ.trieste.it/pinocchio/ 
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Figure 13. Reconstructed PDFs. The upper panels (a,b) show the reconstructed PDFs for the NOG catalogues filled assuming 6rei = 1-2 
(continuous line) compared with the average PDFs computed from the biased NOG catalogues filled with unbiased PSCz simulated 
catalogues (points with errorbars). The left panel (a) shows the case for dmax = 60 Mpc, the right panel (b) shows the case for 
dmax = 30 h^^ Mpc. The lower panels (c,d) show the reconstructed PDFs for the NOG catalogues filled assuming bj-cl = 1 (continuous 
line) compared with the average PDFs computed from the unbiased NOG catalogues filled with unbiased PSCz simulated catalogues 
(points with errorbars). The left panel (c) shows the case for dmax = 60 Mpc, the right panel (d) shows the case for dmax = 30 
Mpc. Dotted lines show a Gaussian with zero mean and unit variance. In each panel the residuals of the reconstructed PDFs with 
respect to that of simulations, rescaled to the variance among the simulations, are shown below. 



clusters, like Virgo, Perseus, Pisces, Centaurus and Hydra, 
are recovered in most if not all cases within a few Mpc. 
We have also checked that the sheet-like overdensity corre- 
sponding to the supergalactic plane is well reproduced. We 
find that, as expected, the NOG reconstruction is better 



than PSCz at high galactic latitude, while it is worse in the 
optical ZOA. 

We conclude that the two catalogues provide comparer 
ble reconstructions of the local Universe. 
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Figure 14. Distribution of the phases of reconstructed initial conditions for the two completions of NOG (continuous lines), compared 
with the initial condition from one simulated catalogue (dotted lines). Fields are padded to zero beyond dmax = 60/1-1 j^jgjj^ 
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5 SUMMARY AND CONCLUSIONS 

In this paper wo have apphed the ZTRACE algorithm 
(Monaco & Efstathiou, 1999) to the NOG (Giuricin et al., 
2000) and PSCz (Saunders et al. 2000) catalogues to recon- 
struct the initial condition of our Local Universe. The use 
of two different catalogues allows us to test and improve the 
reconstruction; in particular, the optical catalogue presents 
the advantages of denser sampling and better relation be- 
tween optical light and galaxy mass, while the infra-red cat- 
alogue benefits from more homogeneous selections and larger 
sky coverage. 

We have studied the relative bias between these cata- 
logues, finding a good linear relation between the density 
fields traced by them, corresponding to a relative bias of 
6rei = 1.1 ± 0.1, in agreement with the expectation based 
on the ratio of the ag parameters for optical and infra-red 
galaxies. The scatter around this relation is found to be 
laxger than that implied by shot noise. This may be a sign 
of stochastic bias; in this hypothesis we obtain a value ~ 0.25 
for the parameter CT6 defined by Dekel & Lahav (1999). 

We have quantified the accuracy of the reconstruction 
through extensive testing with simulated galaxy catalogues. 
The analysis has been performed in the Fourier space. The 
initial conditions reconstructed by ZTRACE reproduce on 
large scales both the fiuctuations of the modules of the 
Fourier modes around the mean (the power spectrum) and 
their phases. The accuracy of the reconstruction drops at 
a scale corresponding to knm — 0.4 h/Mpc, beyond which 
all information is lost. This fenm depends only weakly on 
catalogue, bias scheme, filling of the ZOA and smoothing 
scale. In particular, while for a smoothing radius of ^ 5 — 6 

Mpc the loss of information is determined mainly by 
smoothing, the reconstruction at smaller scales is limited by 



non-linearity. Moreover, the denser sampling of the NOG 
catalogue gives only a modest improvement to the recon- 
struction, which is anyway compensated by the errors in- 
duced by the procedure to fill the ZOA. 

Consistently with the PSCz case (Monaco et al. 2000), 
the initial conditions of the NOG catalogue are found to 
be consistent with the Gaussian, with no convincing dis- 
tortion beyond the flattening of the high peaks induced 
by ZTRACE. We have also checked the distribution of the 
phases, which is found to be consistent with a uniform one 
both for PSCz and NOG. 

The reconstructed initial conditions have been corrected 
for smoothing and peak flattening, using corrections cal- 
ibrated on the simulated catalogues. From these we have 
generated sets of initial conditions suitable for N-body sim- 
ulations by randomly generating the modes beyond feum with 
a power spectrum given by the cosmology in section 3.1. We 
have used the PINOCCHIO tool to run three sets of 10 low- 
resolution simulations, one set from the PSCz and two sets 
for the NOG (with the ZOA filled assuming two different val- 
ues for the relative bias between NOG and PSCz, namely 
6rci = 1 and 1.2). In all cases the large-scale structure ob- 
served in our local Universe is correctly reproduced. The 
locations of the main neighbouring clusters axe also repro- 
duced within a few Mpc. We notice that the NOG-based 
realizations are more accurate than PSCz at high galactic 
latitude, while the opposite is true when the optical ZOA is 
approached. 

With respect to previous works, we improve in many 
regards. We use different galaxy catalogues, NOG and PSCz, 
in place of the smaller 1.2 Jy one (Fisher et al. 1995), used 
by Kolatt et al. (1996) and Mathis et al. (2002); Klypin et al. 
(2001) used the MARK III catalogue of peculiar velocities 
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panels show 3 realizations based on NOG with b^ei = 1.2; mid 
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Figure 16. Observed structures along the supergalactic plane. 
Panel show the location of NOG galaxies. 



(Willick et al., 1997), which allows a reconstruction based 
on linear theory, but at the cost of a much noisier set of 
data. Moreover, we use ZTRACE to self-consistently obtain 
the initial density field from the redshift-space catalogues, 
and then correct for peak flattening and smoothing with 
procedures carefully calibrated on simulations. 
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